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ABSTRACT 

Cesium-based  atomic  standards  are  typically  accurate  to  about 
seven  parts  in  ten  to  the  twelfth.  However,  certain  applications 
require  better  agreement  among  standards.  Qie  exairple  is  that  of  a 
communications  network.  Network-wide  syntonization  requires  each 
standard  in  the  network  to  be  adjusted  to  agree  with  a  master  external 
reference.  This  is  normally  acconplished  by  adjusting  the  "C"  field 
of  each  atonic  standard.  However,  causing  precise  output  frequency 
changes  in  this  way  is  notoriously  difficult.  It  also  reduces  fre¬ 
quency  stability  until  the  standard  eventually  settles  into  a  new 
free-running  frequency. 

There  are  two  ways  in  v^ich  the  output  of  a  cesium  standard  can 
be  adjusted  without  changing  its  "C"  field.  The  internal  synthesizer 
of  the  standard  can  be  adjusted,  or  a  phase  microstepper  can  be  used 
on  its  output.  The  former  technique  shifts  the  output  frequency  in 
small,  discrete  steps,  vhile  the  latter  inserts  small,  discrete  phase 
shifts  into  the  output  signal.  Using  common-view  Global  Positioning 
System  (GPS)  as  a  master  reference,  this  investigation  tests  both 
techniques  and  determines  the  degree  of  syntonization  and  synchroniza¬ 
tion  possible  with  each. 
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INTRODUCTION 


An  important  ancillary  usage  of  the  Global  Positioning  System 
(GPS)  is  that  in  which  it  serves  to  distribute  time  and  frequency  to 
remote  clocks  and  frequency  standards .  In  the  case  in  which  two 
ground-based  clocks  view  the  same  GPS  satellite,  the  theoretical  limit 
on  synchronization  is  about  10  ns  [!]•  GPS  receivers  typically  pro¬ 
vide  synchrcnization  to  GPS  time  better  than  50  ns,  vAiich  neans  that  a 
common-view  GPS  should  allow  synchronization  of  remote  receivers  to 
within  about  100  ns  of  each  other  [1»2]. 

However,  using  the  common-view  ^proach  to  GPS  to  syntonize 
remote  frequency  standards  (i.e.,  to  cause  their  frequencies  to  agree) 
proves  to  be  a  more  challenging  problem.  The  crux  of  the  difficulty 
lies  in  the  relatively  high  effective  frequency  instability  of  the 
corrected  time  transmitted  by  the  satellite.  Figure  1  shows  that 
although  the  long-term  (50  days)  frequency  stability  a  GPS  satellite 
is  about  a  part  in  10^'^,  the  stability  over  one  day  approaches  a  part 
in  10^2.  Suppose  that  we  intend  to  syntonize  two  remote  frequency 
standards  by  measuring  their  difference  to  GPS,  then  adjusting  their 
frequencies  to  reduce  this  difference,  then  repeating  the  process.  If 
the  standards  are  to  be  syntonized  to  within  a  part  in  10^^,  then  each 
iteration  of  the  control  loop  would  require  several  days.  It  would 
take  that  long  to  measure  the  difference  between  the  frequency  stan¬ 
dard  and  GPS  to  the  required  accuracy. 

The  purpose  of  this  note  is  to  determine  vAiether  there  exists  a 
syntonization  method  using  the  commcn-view  GPS  approach  that  gives 
good  results  over  the  short-term  (less  than  1  day) .  The  idea  is  to 
slave  two  frequency  standards  to  the  same  satellite,  then  see  how  well 
the  standards  continue  to  agree  after  the  satellite  is  no  longer  in 
view.  If  identical  control  algorithms  are  used  to  control  identical 
frequency  standards,  using  information  provided  by  identical  GPS 
receivers  viewing  the  same  satellite,  it  seems  reasonable  to  expect 
good  agreement  between  the  frequency  standards. 
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EXPERIMEISITAL  APPARATUS 


The  experiment  was  performed  in  two  phases.  The  experimental 
setup  for  Phase  1  is  shown  in  Figure  2.  In  Phase  1  an  Austron  Model 
2101  GPS  Receiver  measures  the  phase  difference  between  the 
1-pulse-per -second  (1  PPS)  signal  from  the  GPS  satellite  and  that  pro¬ 
duced  by  an  Austron  Model  2310  Disciplined  Cesium  Standard.  This  dif¬ 
ference  is  read  over  the  IEEE-488  bus  by  a  Hewlett  Packard  HP-85 
computer.  A  Kalman  filter  algorithm  running  on  the  HP-85  decides  the 
degree  to  v^ich  the  disciplined  cesium  standard's  frequency  should  be 
shifted  in  order  to  minimize  the  phase  difference.  The  Kalman  filter 
also  strives  to  syntcxiize  the  frequency  standard  to  the  effective  5 
MHz  frequency  oorrespx:9iding  to  the  satellite's  1  PPS  signal.  After 
the  Kalman  filter  calculates  the  appropriate  frequency  shift,  the 
HP-85  commands  (via  the  bus)  the  disciplined  cesium  standard  to  shift 
its  output  frequency  by  this  amount.  The  disciplined  cesium  standard 
is  a  new  device  whose  frequency  can  be  shifted  essentially  instan¬ 
taneously  without  adjusting  the  "C"  magnetic  field  in  its  beam  tube. 

This  setup  was  duplicated  elsewhere  in  the  laboratory.  A  time 
interval  counter  measured  the  phase  difference  between  the  1  PPS  out¬ 
put  of  cesium  standard  in  the  first  setup  and  that  of  the  cesium  stan¬ 
dard  in  the  second.  A  chart  recorder  recorded  this  difference. 

The  setups,  vhich  shall  henceforth  be  referred  to  as  Side  A  and 
Side  B,  shared  the  same  antenna. 

Figure  3  shows  the  experimoital  setup  for  Phase  2.  The  dif¬ 
ference  between  this  setup  and  that  of  Ihase  1  is  that  the  frequency- 
shifting  ability  of  the  disciplined  cesium  standard  is  not  used. 
Instead,  an  Austron  Model  2055  Phase  Microstepper  is  used  to  create  a 
signal  vAiose  frequency  is  shifted  from  that  of  the  input  signal  from 
the  cesium  standard.  The  phase  microstepper  is  a  device  that  shifts 
the  frequency  of  an  input  signal  by  inserting  or  removing  phase  steps 
at  a  rate  spiecified  by  the  user.  The  phase  microstepper  does  not  pro¬ 
vide  a  1  PPS  output.  Hence  the  phase  microstepper's  frequency  output 
was  used  to  drive  a  precision  1  PPS  synthesizer.  As  in  Phase  1  the 
GPS  receiver  informs  the  HP-85  coitputer  of  the  phase  difference 
between  the  output  of  1  PPS  synthesizer  and  the  satellite.  The  Kalman 
filter  calculates  an  amount  by  which  the  frequency  of  the  phase 
microstepper  should  shift  the  frequency  of  the  cesium  standard.  The 
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HP-85  then  commands  the  phase  microstepper  to  produce  this  frequency 
shift.  Because  the  phase  microstepper  requires  BCD  rather  than 
IEEE-488  remote  inpjut,  an  IBM  PC  with  an  IEEE-488  card  was  used  to 
translate  HP-85  IEEE-488  bus  oommands  to  BCD. 

The  Phase  2  setup  was  duplicated  elsewhere  in  the  laboratory.  As 
in  Phase  1,  the  difference  between  the  1  PPS  from  the  two  setups  (Side 
A  and  Side  B)  was  monitored  cn  a  time  interval  counter  and  recorded  cn 
a  chart  recorder. 

The  GPS  receivers  of  Side  A  and  Side  B  shared  the  same  antenna. 

In  both  phases  the  duplication  of  the  HP-85  system  controller  was 
simulated.  That  is,  Side  A  and  Side  B  were  cmtrolled  by  the  same 
computer  that  simultaneously  ran  indep)endent  programs  to  control  each 
side.  In  Phase  2  the  duplication  of  the  IBM  PC  was  simulated.  Side  A 
and  Side  B  shared  the  same  computer  to  translate  IEEE-488  bus  oommands 
to  BCD  instructions  for  their  phase  microsteppers. 

In  both  experimental  phases,  the  receivers  tracked  satellite 
vehicle  13  for  four  hours,  then  stepped  tracking.  Tracking  resumed 
one  hour  later.  The  phase  accumulated  between  Side  A  and  Side  B  was 
recorded.  Tracking  continued  for  thirty  minutes,  then  ceased. 
Tracking  resumed  sixteen  hours  later.  The  accumulaticn  of  phase  bet¬ 
ween  Side  A  and  Side  B  that  occrurred  during  this  sixteen-hour  pericxl 
was  also  recorded.  The  frecjuency  offsets  corresponding  to  the  phase 
accumulation  that  cocurred  during  the  times  that  the  satellite  was  not 
being  tracked  (i.e.,  during  the  "coast"  periods)  were  calculated  and 
are  listed  in  Tables  2  and  4. 


RESULTS 

PHASE  1 

Table  1  sho?s  the  results  of  ten  consecutive  trackings  of 
satellite  13.  While  the  satellites  were  being  tracked  by  the 
receivers,  the  average  of  the  difference  between  the  Side  A  and  GPS  is 
9  ns  with  a  standard  deviation  of  14  ns. 
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Table  2  shows  that  Side  A  and  Side  B  agreed  with  each  other  about 
two  times  itore  closely  than  Side  B  agreed  with  GPS  (the  ratio  of  the 
standard  deviations  was  1.9).  Although  the  difference  between  Side  A 
and  GPS  was  not  monitored  as  closely  as  the  difference  between  Side  B 
and  GPS  (i.e.,  readings  were  not  confirmed  with  a  chart  recorder),  it 
appeared  that  the  difference  between  Side  A  and  GPS  was  ponparable  in 
magnitude  to  the  B-GPS  measurements  listed  in  Table  2.  In  other 
words,  Side  A  and  Side  B  agreed  with  each  other  more  closely  than 
either  agreed  with  GPS. 


Table  1 


MJD 

Phase 

MESSAGE  PGE 

hr 

1  Tracking  Data 

<A-GPS> 

ns 

FMS 

ns 

46761 

6 

3.6 

5 

46761 

5 

25.8 

8 

46760 

6 

2.9 

5 

46760 

6 

39.5 

7 

46759 

2 

3.2 

7 

46759 

6 

-2.3 

4 

46759 

3 

12.2 

8 

46758 

2 

10.7 

62 

46758 

7 

3.2 

5 

46758 

5 

-10.6 

11 

Average 

=  8.8  ns 

Standard  Dev. 

=  1 3 . 7  ns 
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Table  2 


Phase  1  Coast  Data 


A-B 

(frequency  offset) 


5.1  X  10-13 

6.6  X  10-13 

-4.1  X  10-13 

5.7  X  10-13 
-9.6  X  10-13 

9.1  X  10-12 

6.7  X  10-13 

7.8  X  10-12 

7.2  X  10-13 
-9.1  X  10-13 

Avg.  =  1.8  X  10-12 

Standard 

Dev.  =  3 . 4  X  10-12 


B-GPS 

(frequency  offset) 


2.0  X  10-12 
-2.9  X  10-12 

3.2  X  10-12 
-3.9  X  10-12 
-8.1  X  10-12 

7.4  X  10-12 

1.1  X  10-11 
-3.7  X  10-12 

3.4  X  10-12 

1.2  X  10-11 

Avg.  =  2.0  X  10-12 
Standard 

Dev.  =  6 . 4  X  10-12 


PHASE  2 

Table  3  shows  that  the  phase  microstepper  was  able  to  track  the 
GPS  an  order  of  magnitude  more  closely  than  the  disciplined  cesium 
standard  used  in  Phase  1.  Ihe  average  phase  difference  was  0.8  ns 
with  a  stcindard  deviaticxi  of  1.5  ns,  although  the  sairple  size  is 
small.  This  iirproved  synchronization  may  be  due  to  microstepper's 
ability  to  cause  frequency  shifts  ten  times  larger  than  the  maximum 
frequency  shifts  that  can  be  produced  by  the  disciplined  cesium  stan¬ 
dard. 


At  first  glance  of  Table  4  it  ^p)ears  as  though  syntcnization  in 
Phase  2  was  also  superior  by  an  order  of  magnitude:  -2.3  x  10-13 
instead  of  1.8  x  10-12.  However,  the  standard  deviation  of  2  x  10-12 
is  much  larger  than  the  sanple  mean  2.3  x  10-13.  Thus  the  relevant 
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Table  3 


Phase  2  Tracking  Data 


WJD 

MESSAGE  AGE 

hr 

<A-GPS> 

ns 

PMS 

ns 

46764 

6 

0.8 

5 

46764 

4 

2.7 

57 

46763 

6 

-0.1 

5 

46763 

6 

-0.6 

4 

46763 

3 

3.5 

58 

46762 

6 

-0.8 

4 

46762 

6 

.3 

7 

Average  =  0.8  ns 
Standard  Dev.  =  1.5  ns 


Table  4 

Fbase  2  Coast  Rate 
A-B  B-GPS 


_ (frequency  offset) 

4.2  X  10-13 
-7.3  X  10-13 
-4.6  X  10-12 

8.0  X  10“13 

7.2  X  10-13 
-6.1  X  10-13 

2.4  X  10-12 

Avg.  =  -2.3  X  10“13 
Standard 

Dev.  =  2.0  X  10-12 


( frequency  offset) 

-1.0  X  10-12 

-2.9  X  10-12 

1.2  X  10-11 

-3.9  X  10"12 
-6.2  X  10-12 
9.4  X  10-12 

1.2  X  10-11 

Avg.  =  2.8  X  10-12 

Standard 

Dev.  =  7.4  X  10-12 
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measure  of  syntonization  in  this  case  is  the  standard  deviation,  not 
the  average  offset.  Hence,  the  syntonization  achieved  in  Phase  2  (2.0 
X  ig  conparable  to  that  achieved  in  Phase  1  (3.4  x 

As  in  Phase  1,  Side  A  and  Side  B  were  syntonized  to  each  other 
(standard  deviation  of  2.0  x  10~12)  tetter  than  they  were  syntonized 
to  GPS  (standard  deviation  of  7.4  x 

aaJCLUSiasFS 

While  a  satellite  is  in  common  view,  it  is  possible  to  use  a 
Kalman  filter  to  synchronize  disciplined  cesium  standards  to  each 
other  to  better  than  15  ns.  When  the  frequency  shifting  is  performed 
ty  phase  microsteppers  instead  of  by  the  disciplined  cesium  standards, 
syntonization  was  about  an  order  of  magnitude  closer.  This  irey  be  due 
to  the  microstepper's  ability  to  remove  large  phase  differences  (which 
were  usually  dbserved  at  the  end  of  the  interval  during  which  the 
satellite  was  not  being  tracked)  at  a  rate  ten  times  faster  than  the 
disciplined  cesium  standards. 

The  syntonizations  achieved  between  disciplined  cesium  standards 
was  about  4  x  10“^ 2.  However,  the  syntonization  between  the  cesium 
standards  was  twice  as  close  v^en  the  frequency  shifting  was  performed 
by  fhase  microsteppers.  Again  this  may  be  due  to  the  microstepper's 
ability  to  itove  more  rapidly  than  the  disciplined  cesium  standard.  In 
both  cases,  when  the  satellite  was  not  in  view,  the  syntonization 
between  the  independent  frequency  standards  to  each  other  was  better 
(by  a  factor  of  2  to  4)  than  the  syntonization  of  the  standards  to 
GPS.  This  is  probably  due  to  the  relatively  poor  short-term  stability 
of  GPS  caused  by  upload  transmissicns. 

Although  in  the  current  C3^S  environment  the  phase  microstepper 
performs  somewhat  better  than  the  disciplined  cesium  standard  in 
achieving  short-term  syntonizaticxis  and  synchrcnizations,  it  seems 
likely  that  the  disciplined  cesium  standard  is  better  suited  to  the 
long  term.  This  is  because  the  disciplined  cesium  standard  has  a 
theoretical  frequency-shift  resolution  70  percent  better  than  that  of 
the  phase  microstepper.  As  GPS  is  iirproved  and  its  effective  fre¬ 
quency  instability  is  reduced,  sole  advantage  of  the  phase 
microstepper  approach  (faster  slew  rate)  will  be  cutweighted  by  its 
disadvantages  —  the  need  for  additicnal  hardware  and  the  reduction  of 
spectral  purity. 
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In  summary,  the  disciplined  cesium  standard  provides  an  alter¬ 
native  means  of  syntonizing  to  a  master  reference  (such  as  GPS).  Its 
effectiveness  relative  to  the  phase  microstepper  approach  is  limited 
only  by  its  slower  slew  rate.  Ihis  makes  it  well  suited  for  Icng-terra 
applications  involving  extremely  stable  master  references. 
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USNO  data  reprinted  from  Advanced  GPS  Receiver  Design  by 
Navigation  Technology  Seminars,  Inc.,  Falls  Church  1984 

Figure  1. 
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INTERVAL 

COUNTER 


In  actual  setup  Side  A  and 


ide  B  shared  HP-85,  IBM-PC,’  and  GPS  antenna. 


Figure  3 


QUESTIONS  AND  ANSWERS 


JOHN  DICK,  JET  PROHJLSION  LABORATORY:  What  you  would  expect  from  this  sort  of 
system  depends  on  the  algorithm  that  you  have  in  the  computer.  Can  you  charac¬ 
terize  that  in  terms  of  an  integration  time  or  something  else,  besides  the 
number  which  you  did  give,  the  slew  rate. 

MR.  RODRIGO:  It  is  a  two  state  Kalman  filter  algorithm  designed  by  Jim  Barnes. 
The  time  constant  parameters  in  it  are  taken  to  be  very  large,  corresponding 
to  short  time  constants  at  the  start.  The  algorithm  just  goes  about  learning 
what  an  appropriate  time  constant  is.  When  it  starts,  you  are  talking  about 
time  constants  on  the  orxier  of  the  data  interval,  15  seconds  or  so.  By  the 
time  it  is  finished  the  algorithm  has  decided  what  an  appropriate  time  con¬ 
stant  is.  I  don’t  really  have  an  answer  for  v4iat  that  number  is. 


